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Abstract. Non-LTE line formation calculations of Fei are performed for a small number of reference stars to investigate and 
quantify the efficiency of neutral hydrogen collisions. Using the atomic model that was described in previous publications, the 
final discrimination with respect to hydrogen collisions is based on the condition that the surface gravities as determined by the 
Fe i/Fe n ionization equilibria are in agreement with their astrometric counterparts obtained from Hipparcos parallaxes. 
High signal-to-noise, high-resolution echelle spectra are analysed to determine individual profile fits and differential abun- 
dances of iron lines. Depending on the choice of the hydrogen collision scaling factor Sh, we find deviations from LTE in Fei 
ranging from 0.00 (Sh = °°) to 0.46 dex (Sh = for HD140283) in the logarithmic abundances while Fe n follows LTE. With the 
exception of Procyon, for which a mild temperature correction is needed to fulfil the ionization balance, excellent consistency is 
obtained for the metal-poor reference stars if Balmer profile temperatures are combined with Sh = 3. This value is much higher 
than what is found for simple atoms like Li or Ca, both from laboratory measurements and inference of stellar analyses. 
The correct choice of collisional damping parameters ("van-der-Waals" constants) is found to be generally more important for 
these little evolved metal-poor stars than considering departures from LTE. For the Sun the calibrated value for Sh leads to 
average Fe i non-LTE corrections of 0.02 dex and a mean abundance from Fe i lines of log e(Fe) = 7.49 + 0.08. 
We confront the deduced stellar parameters with comparable spectroscopic analyses by other authors which also rely on the 
iron ionization equilibrium as a gravity indicator. On the basis of the Hipparcos astrometry our results are shown to be an order 
of magnitude more precise than published data sets, both in terms of offset and star-to-star scatter. 

Key words. Line: formation - Sun: abundances - Stars: fundamental parameters - Stars: abundances - Stars: individual: HD 
140283 - Stars: individual: Procyon 



1. Introduction 

As is evident from our previous investigations (Gehren et al. 
2001a, 2001b), the large scatter among laboratory/ values of 
neutral iron lines does not allow the unambiguous determina- 
tion of a kinetic model of the iron atom in the solar atmosphere. 
In particular, the role of collisions with neutral hydrogen atoms 
could not be fully explored. This resulted in an uncertainty of 
the hydrogen collision scaling factor Sh (a multiplicative mod- 
ification to the Drawin (1968, 1969) formula for allowed in- 
elastic collisions with neutral hydrogen) which becomes im- 
portant as soon as the kinetic model is applied to metal-poor 
stars. Since in the atmospheres of halo and thick disk stars 
line-blanketing in the UV becomes more and more unimportant 
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leading to an increase of photoionization, while simultaneously 
electron collisions can no longer counterbalance the predomi- 
nance of radiation interactions, the hydrogen collisions are an 
essential part of the kinetic equilibrium of iron, and many other 
photoionization-dominated metals as well. This problem is fur- 
ther enhanced by the new radiative cross-sections of Bautista 
(1997) which are substantially larger than simple hydrogenic 
approximations would make believe. 

To establish the influence of deviations from LTE on the 
ionization equilibrium of iron we have therefore chosen a num- 
ber of reference stars which we feel are sufficiently represen- 
tative of the local metal-poor stellar generation on and near the 
main sequence. Our interest is particularly devoted to the im- 
portant question: to what extent is the derivation of spectro- 
scopic stellar parameters and metallicities affected by the ki- 
netic equilibrium of iron? Answering this question involves the 
investigation of the radiation field (depending mostly on the ef- 
fective temperature) and of the collisional interaction (depend- 
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Table 1. Basic stellar parameters of the reference stars. The iron abun- 
dance is derived from Fe n lines with the microturbulence fulfilling the 
usual line strength constraint (see text). [a/Fe] = 0.4 was adopted in 
the computation of the metal-poor atmospheres, for Procyon (= HD 
61421) a solar mixture was assumed. Except for Procyon (data from 
Allende Prieto et al. 2002), surface gravities are from Hipparcos paral- 
laxes with masses derived from the tracks of VandenBerg et al. (2000). 
[Fe/H] = log(Fe/H), - log(Fe/H) 







logg 


[Fe/H] 




tthip 


M 


star 


[K] 


(cgs) 




[km/s] 


["] 


[M Q ] 


HD 103095 


5070 


4.69 


-1.35 


0.95 


109.2 


0.62 


HD 19445 


6032 


4.40 


-2.08 


1.75 


25.9 


0.67 


HD 140283 1 


5806 


3.69 


-2.42 


1.70 


17.4 


0.79 


HD 84937 1 


6346 


4.00 


-2.16 


1.80 


12.4 


0.79 


Procyon 


6510 


3.96 


-0.03 


1.83 


285.9 


1.42 


Sun 


5780 


4.44 


0.0 


1.0 




1.00 



1 for this star extinction due to interstellar reddening was considered 
following Hauck & Mermilliod (1998, cf. Table 7) 



ing mostly on surface gravity). Therefore our choice of metal- 
poor reference stars has led to four representatives: HD 103095 
(= Gmb 1830) as a cool main sequence star, HD 19445 as a 
typical subdwarf, HD 84937 as a turnoff star and HD 140283 
as a moderately cool subgiant. The last two of these reference 
stars have also been chosen by Asplund et al. (2000) to model 
hydrodynamic convection in metal-poor stars. 

Photospheric surface gravities are usually determined by 
means of three independent methods. In principle the most re- 
liable approach uses astrometry (presently the Hipparcos par- 
allaxes, ESA 1997) with an appropriate guess of the stellar 
mass, which enters the resulting logg only with a cr(logM) 
uncertainty. Since the mass uncertainty of these very old stars 
is generally below 0.05 M (when referring to a given set of 
evolutionary tracks), the corresponding error stays below 0.03. 
The uncertainty of the effective temperature enters as 4 <x(log 
T e ff) ~ 0.02. Thus the prime error in log g is due to that of the 
parallax itself (2 <x(log n)), which for the four stars is between 
0.01 (Gmb 1830) and 0.07 (HD 84937). The maximum error 
then is between 0.06 and 0.12 in logg. 

The second method refers to the damping wings of strong 
neutral metal lines such as the Mgib triplet. However, as 
demonstrated by Fuhrmann (1998a), for stars with [Fe/H] < -2 
the wings of such lines become extremely shallow and the 
determination correspondingly uncertain. The combination of 
Balmer profile temperatures and Mgib gravities (both in LTE) 
was shown to be fully compatible with the Hipparcos astrome- 
try (Fuhrmann 1998b, 2003). 

Originally, the surface gravity of remote subdwarfs and 
subgiants was derived from the ionization equilibrium of iron 
where for simplicity and lack of better evidence one used 
the assumption that Fei/Fen is populated according to the 
Saha/Boltzmann equilibrium (LTE). For a number of stars in- 
cluding Procyon as a primary standard this was shown to lead 
to results incompatible with those obtained from strong line 
wings (Fuhrmann et al. 1997) or Hipparcos parallaxes. More 
specifically, evolutionary stages are predicted which have not 



Table 2. Comparison of effective temperature determinations for the 
five reference stars. TI99 = Thevenin & Idiart (1999), F00 = Fulbright 
(2000), AAM-R96 = Alonso et al. (1996), GCC96 = Gratton et al. 
(1996), KSG03 = this work. The mean offset is computed for the four 
metal-poor stars according to A T c ff (study - KSG03). Significant off- 
sets with respect to our temperature scale are encountered in the stud- 
ies of Thevenin & Idiart (1999) and Fulbright (2000) 
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yet been reached. Therefore, the present approach has aban- 
doned the LTE assumption and obtains the surface gravity from 
a fit to the Fe lines in kinetic equilibrium. 

To evaluate the reliability of the non-LTE calculations in 
stars with different parameters Procyon (s a CMi) was also 
selected, because the discrepancy between logg obtained from 
ionization equilibrium and parallax remains unresolved. 

Table 1 lists the basic stellar parameters of the five stars 
plus the Sun. 

2. Observations 

Spectra of the stars were taken with the fibre-fed echelle spec- 
trograph FOCES (Pfeiffer et al. 1998) at the 2.2m telescope of 
the German-Spanish Astronomical Centre on Calar Alto dur- 
ing a number of observing runs between 1999 and 2001, with 
a spectral resolution of R = 65 000 (HD 140283 was observed 
at R = 40 000). As measured in line-free spectral windows near 
Ha, the peak signal-to-noise ratio (S/N) exceeds 300:1 in all 
cases. With a spectral coverage of 4200 - 9000 A our analyses 
naturally focus on the lines of Fe i and Fe n in the visible disre- 
garding the resonance lines in the (near) UV. 
Owing to the design of FOCES (the light of both star and flat- 
field lamp being transmitted through the same optical path), 
the resulting flatfield calibration of the stellar spectra produces 
an extremely predictable continuum run that contains only the 
remaining second-order variation with wavelength which re- 
sults from the difference between stellar T e g and the tempera- 
ture of the halogen lamp. The internal accuracy of the effective 
temperature determination is therefore very high, with profile 
residuals well below 0.5 % of the continuum (cf. Korn 2002). 
The profile analysis of metal lines also profits from this quality, 
especially when the line density is high. 

3. Stellar effective temperatures and gravities 

The temperature scale adopted is based on Balmer line- 
broadening theory combined from Stark effect profiles of Vidal 
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et al. (1970) and the Ali & Griem (1965, 1966) approximation 
of resonance broadening (see Fuhrmann et al. 1993). The basis 
of this approach was the different relative contributions of these 
broadening processes to Ha on the one hand (essentially broad- 
ened by H collisions) and the higher series members (domi- 
nated by the Stark effect) on the other. The results have been 
fixed with a minimum number of free parameters for solar-type 
stars with ranging metallicities based on empirical evidence. 
Weights are assigned to the temperatures derived from each 
line based on temperature sensitivity, S/N and blaze correc- 
tion reliability. At present we refrain from adopting the recent 
hydrogen line-broadening theory of Barklem et al. (2000). As 
our own implementation shows, it does not reproduce the solar 
temperature to within 100K (using the MARCS code and a^ 2 
fitting technique which allows for the intersection of observed 
profiles this value is 50 K, Barklem et al. (2000). In addition, it 
fails to reproduce the profile shape of higher Balmer lines (H/3 
and up) of stars like HD 140283 and HD 19445, both of which 
are located in the temperature band where differences between 
the two broadening theories are expected to be largest. 

We utilize the atmospheric model MAFAGS (Fuhrmann et 
al. 1997) which is line-blanketed by means of the Kurucz ODFs 
(Kurucz 1992). Note that these ODFs were computed assum- 
ing 1 log e(Fe) = 7.67. We therefore rescale them by -0.16dex 
to account for the low Solar iron abundance of log e(Fe) = 7.51. 
Unlike in other widely-used atmospheric codes like ATLAS or 
MARCS, we use a mixing length a = l/H p = 0.5 in the frame- 
work of the Bohm-Vitense (1958) theory of convection in or- 
der to bring Ha and the higher Balmer lines into concordance. 
The differences with respect to other methods of temperature 
determination have been described in Fuhrmann et al. (1993) 
and Fuhrmann (1998b). We note that the choice of the mixing- 
length affects the Balmer lines, but the influence on the Fe i 
and Fe n lines analysed in the Sun turns out to be equivalent 
to A log e Fe < 0.01 dex. The same is true for metal-poor stellar 
atmospheres such as that of HD 140283. It also holds for the 
corresponding non-LTE computations. 

The temperature scale resulting from our hydrogen line 
profile fitting is higher for the metal-poor stars than the photo- 
metric one of Thevenin & Idiart (1999) and Fulbright (2000) 
by more than 100K. It is, however, in excellent agreement 
with the infrared flux method (IRFM) temperatures of Alonso 
et al. (1996) and Gratton et al. (1996). Table 2 confronts the 
corresponding numbers in detail. Only HD 140283 displays a 
markedly lower temperature in the Alonso et al. calibration. 
This is, at least in part, due to the zero-reddening assumption 
adopted by these authors. A discussion of this point can be 
found in Fuhrmann (1998a). 

Yet another independent temperature scale seems to sup- 
port the Balmer profile temperatures: Peterson et al. (2001) 
derive temperatures for HD 19445 and HD 84937 (and other 
stars) from absolute fluxes in the mid-UV. Effective temper- 
atures of 6050 K and 6300 K were obtained. Peterson (2003, 
priv. comm.) analysed HD 140283 with similar methods ob- 
taining 5850 K. In total the mean offset between our tempera- 
ture scales is thus well below 10 K. 

1 log£(Fe) = log(n Fe /n H ) + 12 



For the calibration the surface gravity is taken from 
Hipparcos parallaxes. Additional quantities entering the equa- 
tion are taken from Alonso et al. (1996, bolometric correction 
BC), Hauck & Mermilliod (1998, interstellar extinction A v ) 
and VandenBerg et al. (2000, stellar mass M). To derive the 
latter, the star is placed among the evolutionary tracks appro- 
priately chosen with respect to metallicity ( [Fe/H] nlte) and a- 
enhancement ([0/Fe] NLTE , according to Reetz 1999). This pro- 
cess is iterated to convergence. 

After the calibration the final atomic model (which best re- 
produces the astrometric results on average) is used to derive 
gravities which can then be compared with Hipparcos to quan- 
tify the star-to-star scatter. 

4. Metal-poor reference stars and Procyon 

As described in Papers I and II in considerable detail, our 
present approach to the analysis of the reference stars is based 
on a number of assumptions. 

Atmospheric models are calculated for all stars (includ- 
ing the Sun) according to the same type of input parameters, 
boundary conditions, and constraint equations. We thus inten- 
tionally avoid the use of empirical models. 

The major orientation of our analysis is that of empirical 
evidence. Although some theoretical model descriptions appear 
more satisfactory than others, we build upon the experience that 
shortcomings in theories can be detected only by comparison 
with observations. 

Our basic approach to stellar spectroscopy is differential 
in the sense that we require all types of input physics to be 
in agreement with solar observations. This requirement clearly 
dictates the individual line fit approach, and it excludes the use 
of absolute / values. It is important to recognize that the typ- 
ical abundance scatter obtained from solar spectral analysis of 
Fe lines is considerably higher than warranted by both obser- 
vational error and differences in theoretical models. There are 
essentially two explanations for this result: either the / values 
are not reliable or a significant fraction of Fe lines are contam- 
inated by relatively strong (yet undetected) blends. 

Differential spectrum analysis must account for line broad- 
ening, particularly when solar-type and extremely metal-poor 
stars are to be compared. This is in fact the most impor- 
tant constraint because it involves additional atomic parame- 
ters and line broadening theories. Every strategy to deal with 
strong line wings automatically involves the problems of at- 
mospheric modelling. Consequently, for interpretation we es- 
tablish steps of decreasing reliability with constraints coming 
from the observed line profile and the concept that the solar 
photospheric abundance of Fe is the same for all lines. Purely 
radiative parameters (/ values) are complemented by colli- 
sional broadening (Stark effect, van der Waals-type damping). 
The variety of line strengths observed in the solar spectrum 
then provides a comparison between weak and strong lines of 
the resulting abundances, with all / values taken from labora- 
tory measurements and the damping constants adjusted for the 
strong lines to minimize trends with line strength (see Paper 
II). Although this interferes occasionally with thermal or mi- 
croturbulent broadening velocities, Doppler broadening can be 
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generally separated from the influence of damping wings. As 
a result, stellar line abundances are compared with the individ- 
ual abundances of their solar counterparts. This is what is often 
called a line-by-line differential (non-)LTE analysis. 

The line broadening problems are less significant for a dif- 
ferential analysis of Fe II lines, however, here Doppler broad- 
ening plays an important role. The scatter introduced to the so- 
lar Fe n abundance by the various data sets is even higher than 
for Fe i thus a differential approach is absolutely necessary. 

The above considerations of the different contributions to 
line formation severely limit the possible conclusions. Still, 
experience with stellar analyses show that atomic data are in 
many cases better determined from the solar spectrum than 
from a terrestrial laboratory. 

4.1. Non-LTE calculations 

With little progress in the understanding of the physics of hy- 
drogen collisions, the review of Lambert (1993) still reflects 
our current state of knowledge. While the energy range en- 
countered in stellar atmospheres is not accessible to exper- 
iment, Lambert concludes from measurements at 15eV that 
the Drawin (1968, 1969) formula most likely overestimates 
the true efficiency of this thermalizing process by two orders 
of magnitude. If this was the case, they could safely be dis- 
regarded. Rather than assuming that hydrogen collisions are 
inefficient all together, we take on the view of Steenbock & 
Holweger (1984) that the "Jack of reliable collision cross- 
sections constitutes a major problem in statistical-equilibrium 
studies of cool stars" '. To face this problem we keep the ef- 
ficiency of hydrogen collisions Sh as a free parameter which 
must be determined by appropriate observations. 

As would be predicted, electron collisions cannot con- 
tribute very much to thermalization, since the density of free 
electrons is roughly a factor 10 smaller than in comparable 
stars of solar metal content. Our calculations support this idea 
fully. In metal-poor stars therefore no other process but hy- 
drogen collisions can - at least to some extent and depending 
on their cross-section - compensate for the high photoioniza- 
tion rates which would otherwise lead to extreme overioniza- 
tion at variance with the spectroscopic evidence given below. 
It is of course unsatisfactory to use such a rather primitive rep- 
resentation of the collision process as the Drawin approxima- 
tion which was never meant to be extended to anything else 
but collisions between rare gases. There is, however, a similar 
behaviour of both formulae of van Regemorter (1962, for elec- 
tron collisions) and Drawin (for hydrogen collisions) in that the 
rates for both types of thermalization are proportional to the / 
value and both depend on the inverse threshold energy. Such 
a variation leads to a general increase of collision rates with 
excitation energy, because highly excited levels cluster near 
the ionization limit. Solar photospheric emission lines show at 
least for Mg i and Ah that it is just these highly excited levels 
which lead to Rydberg transitions dominated by a population 
inversion. This requires low collision rates (whether electron 
or hydrogen) for such transitions, and in fact the corresponding 
hydrogen collision factors for such atomic models were found 
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Fig. 1. Departure coefficients b t for Fe i terms in HD 140283 as a func- 
tion of the optical depth log t 50 oo- It becomes evident that they depend 
drastically on Sh- The line source functions are Planckian for all mod- 
els with Sh > 0, but not for the Sh = model. The best model (Sh = 3) 
already returns abundances quite close to LTE 



to fit the observed emission lines with Sh values significantly 
below unity (Baumuller & Gehren 1996, Zhao et al. 1998). 
Carlsson et al. (1992) reproduce the IR emission lines entirely 
without hydrogen collisions, but do not simultaneously model 
the optical absorption lines. 

For Fe, we determine Sh by requiring that both Fe i and 
Fe ii lines produce a single-valued mean iron abundance when 
the gravity is assumed to be that derived from the Hipparcos 
parallax. This leaves the hydrogen collision enhancement fac- 
tor Sh as the single free parameter in this investigation. 

A number of non-LTE test calculations with various 
choices for Sh was performed for the reference stars, together 
with LTE calculations. Line profiles were fitted using exter- 
nal Gaussian broadening profiles essentially defined by the in- 
strumental resolution. While the intrinsic profiles are better 
reproduced by a radial-tangential approximation to macrotur- 
bulence (Gray 1977), the dominance of the instrumental res- 
olution and the relative weakness of the iron lines considered 
(Wi < 100 mA) justify the use of a Gaussian. Equivalent widths 
(EWs) are then computed from the theoretical spectra to con- 
strain the microturbulence Thus EWs only enter the analysis 
for the derivation of all other aspects are performed using 
line profiles. The corresponding data is presented in Tables 4 
and 5. 

4.2. Evidence for non-vanishing hydrogen collisions: 
the subgiant HD 140283 

The dynamic range of Fei abundances in HD 140283 upon 
varying Sh is exemplified by Figure 1 and Table 3: the extreme 
underpopulation of Fe i in its atmosphere is clearly seen. The 
depopulation with respect to the LTE case already starts in deep 
photospheric layers. For vanishing Sh the departure coefficients 
bi drop to 0.5 near log t = -0.5. 
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Table 3. Fei and Fen abundances for the reference stars and the Sun for different assumptions concerning the efficiency of hydrogen collisions. 
Differential abundances [Fe/H]<iifr are given, except in the case of the Sun where log e(Fe) is specified, a refers to the line-to-line scatter which is 
a factor yfn larger than the error of the mean. Note that the scatter in [Fen/H] diff is significantly reduced. Like in the case of HD 140283 discussed 
above, a strong imbalance between Fe i and Fe n results for all metal-poor stars if Sh = is assumed. Instead, the ionization equilibrium points 
towards S H = 3. Procyon cannot be fitted by any model 
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Fig. 2. A comparison between the stellar parameters of HD 140283 
for different assumptions concerning the efficiency of hydrogen col- 
lisions. Sh = 3 corresponds to our best estimate and places this star 
in a subgiant phase of evolution. A mass close to 0.8 M Q is indicated 
from a track of appropriate metallicity by VandenBerg et al. (2000). 
However, assuming Sh = (no collisions) results in a r c ff-log g combi- 
nation in irreconcilable conflict with the star's known properties, both 
in terms of its evolutionary stage and mass 



Models with Sh > produce less extreme departures as ra- 
diative processes (photoionization) are increasingly compen- 
sated for by collisional ones (hydrogen collisions). More im- 
portantly, though, a non-vanishing Sh secures relative thermal 
excitation equilibrium: although strongly overionized with re- 
spect to Fe ii, the Fe i populations are kept thermal relative to 
one another, even by collisional efficiencies lower than indi- 
cated by the Drawin formula (e.g. Sh = 0.1). Since only the 
ratio of the departure coefficients enters the line source func- 
tion S(,, identical bj result in a Planckian S l y . This means that 
within a given model non-LTE corrections only depend on the 
contribution function, i.e. the depth of formation. Thus low- 
excitation lines will experience larger corrections on average. 



Calculations including hydrogen collisions then consti- 
tute a model sequence with non-LTE abundance corrections 
steadily decreasing with an increasing value for Sh- The be- 
haviour found in the metal-poor stars is morphologically simi- 
lar to the case of the Sun. We thus refer the reader to Papers I 
& II for more details. 

Table 3 shows how these departure coefficients translate 
into Fe i abundances. With a lcr scatter of 0.05 dex in Fe n the 
S H = model is ruled out at more than 3cr. In addition to the 
discrepant mean values, the Sh = model would require a sig- 
nificantly smaller microturbulence for Fei (£ Fe i = 1.2km/s), in 
conflict with the one derived from Fe n lines (£ FeU = 1.7km/s). 
If the stellar parameters were determined assuming Sh = (this 
time using Ha as the temperature indicator which shows some 
dependence on gravity), the ionization equilibrium would be 
reached at T eB = 5550K, log g = 4.25 and [Fe/H] = -2.26 (cf. 
Fig. 2). Assuming a sensible stellar mass of M = 0.8M©, A H ip 
(as defined in the caption of Table 7) would be -54 %. 

4.3. The other metal-poor reference stars 

On the Balmer profile temperature scale Gmb 1830, HD 19445 
and HD 84937 all require an Sh close to 3. This is surprising 
as these stars are in different evolutionary phases and require 
different non-LTE corrections. In fact, the non-LTE effects in 
Gmb 1830 are smaller than in the Sun resulting in an abun- 
dance and gravity correction with an opposite sign from all the 
other stars. This turns out to reduce the star-to-star scatter sig- 
nificantly (cf. Fig. 4). Non-LTE effects for HD 19445 are as 
large as in the Sun. The differential non-LTE result thus coin- 
cides with the LTE one. 

With S H set to 3 the non-LTE corrections are generally 
quite small and never entail a gravity correction larger than 
0. 1 dex. It might come as a surprise that we are nonetheless 
able to remove the large discrepancies between Mgib and 
Fe i/Fe n as gravity indicators uncovered by Fuhrmann (1998a). 
He points out that using Fei/Fen results in gravities up to 
0.5 dex smaller than those from Mgib which are in excellent 
agreement with Hipparcos. The solution to this apparent con- 
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Table 4. Atomic data of Fe i lines used in the analysis of the reference stars. Column 1 gives the multiplet number, column 2 the wavelength 
in A, column 3 the lower excitation energy, column 4 the differential logg/ value with respect to log£(Fe) = 7.51, column 5 the employed 
log C 6 value. The products log gfs Q that the log gf values presented here imply differ slightly from those presented in Paper II. This is because 
a global minimization was performed in Paper II to determine the logC6 values, whereas lines were fitted individually here on the basis 
of the minimization procedure. A radial-tangential approximation to macroturbulence (Gray 1977) was used for the analysis of the Sun, a 
Gaussian otherwise. The equivalent widths (in mA) in columns 6-10 are of low quality, but only enter the analysis in the derivation of the 
microturbulence. We discourage from using these values for abundance analyses 

mult. A E ]ow logg/ d i ff logCg Procyon Gmb 1830 HD 19445 HD 84937 HD 140283 



1 


5166.2 


0.000 


-4.18 


-32.22 


73 




16 


5 


1 1 


1 


5225.5 


0.110 


-4.74 


-32.20 


33 


62 








1 


5247.0 


0.087 


-4.91 


-32.21 


27 


57 








1 


5250.2 


0.121 


-4.88 


-32.20 


26 


57 








2 


4347.2 


0.000 


-5.47 


-32.16 




27 








2 


4427.3 


0.052 


-2.94 


-32.15 






60 


36 


56 


2 


4445.4 


0.087 


-5.43 


-32.15 


11 


27 








3 


4232.7 


0.110 


-4.96 


-32.12 




47 








13 


6498.9 


0.958 


-4.64 


-32.08 




27 








15 


5269.5 


0.859 


-1.34 


-32.04 






111 


84 


97 


15 


5328.0 


0.915 


-1.49 


-32.03 






100 


74 


89 


15 


5371.4 


0.958 


-1.70 


-32.02 






89 


64 


80 


15 


5397.1 


0.915 


-2.03 


-32.03 






73 


49 


64 


15 


5405.7 


0.990 


-1.84 


-32.02 






76 


51 


69 


34 


6581.2 


1.485 


-4.73 


-31.96 




7 








34 


6739.5 


1.557 


-4.92 


-31.91 




4 








36 


5194.9 


1.557 


-2.15 


-31.83 


109 




34 


17 




36 


5216.2 


1.608 


-2.22 


-31.82 


101 




27 


13 


22 


41 


4404.7 


1.557 


-0.18 


-31.71 








98 




42 


4147.6 


1.485 


-2.15 


-31.67 


100 






17 


26 


42 


4271.7 


1.485 


-0.25 


-31.70 








103 




62 


6151.6 


2.176 


-3.31 


-31.72 


20 


22 








62 


6297.8 


2.223 


-2.70 


-31.72 


46 


49 








64 


6082.7 


2.223 


-3.60 


-31.70 




13 








64 


6240.6 


2.223 


-3.30 


-31.71 


21 


21 








66 


5079.2 


2.198 


-2.01 


-31.58 


92 


93 








66 


5198.7 


2.223 


-2.14 


-31.59 


77 










66 


5250.6 


2.198 


-2.00 


-31.61 


90 


87 








68 


4494.5 


2.198 


-1.17 


-31.45 










36 


69 


4442.8 


2.176 


-2.79 


-31.44 




37 








71 


4282.4 


2.176 


-1.06 


-31.39 






57 


41 


45 


111 


6421.3 


2.279 


-1.98 


-31.71 


87 


91 








111 


6663.4 


2.424 


-2.42 


-31.67 


60 


54 








111 


6750.1 


2.424 


-2.56 


-31.68 


45 


46 








111 


6978.8 


2.484 


-2.44 


-31.67 


53 










114 


4924.7 


2.279 


-2.15 


—31 59 


73 










114 


5049.8 


2.279 


-1.31 


-31.54 






33 


18 


25 


114 


5141.7 


2.424 


-2.22 


-31.50 


62 










116 


4439.8 


2.279 


-3.04 


-31.39 


24 


24 








152 


4187.0 


2.449 


-0.66 


-30.79 






62 


43 


54 


152 


4222.2 


2.449 


-0.98 


-30.80 






42 


26 


34 


152 


4233.6 


2.482 


-0.70 


-30.79 






59 


40 


44 


152 


4250.1 


2.469 


-0.49 


-30.81 






72 


48 


53 


152 


4260.4 


2.399 


-0.01 


-30.84 






100 


74 


75 


168 


6393.6 


2.433 


-1.46 


-31.65 


101 




22 


13 


15 


168 


6494.9 


2.404 


-1.25 


-31.67 






38 


20 


28 


168 


6593.8 


2.433 


-2.32 


-31.66 


54 


57 








169 


6136.6 


2.453 


-1.39 


-31.62 






26 


15 


18 


169 


6191.5 


2.433 


-1.46 


-31.63 






25 


13 


18 


169 


6252.5 


2.404 


-1.62 


-31.64 


95 




18 


11 


12 


170 


5916.2 


2.453 


-2.90 


-31.59 


25 
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continuation of Table 4 



mult. 


A 


^low 


logg/diff 


logC 6 


Procyon 


Gmb 1830 


HD 19445 


HD 84937 


HD 140283 


206 


6609.1 


2.559 


-2.62 


-31.62 


34 


33 








207 


6065.4 


2.608 


-1.46 


-31.55 


97 


99 


16 


7 


11 


207 


6137.6 


2.588 


-1.31 


-31.57 










15 


207 


6200.3 


2.608 


-2.35 


-31.57 


44 


41 








207 


6230.7 


2.559 


-1.16 


-31.59 






30 


15 


21 


207 


6322.6 


2.588 


-2.35 


-31.58 


48 


44 








209 


5778.4 


2.588 


-3.51 


-31.53 


7 










268 


6546.2 


2.758 


-1.59 


-31.54 


81 


88 






8 


268 


6592.9 


2.727 


-1.47 


-31.55 




92 






9 


268 


6677.9 


2.692 


-1.33 


-31.57 


100 




20 


11 


12 


318 


4890.7 


2.875 


-0.44 


-30.80 






50 


33 


38 


318 


4891.4 


2.851 


-0.17 


-30.81 






67 


48 


54 


318 


4918.9 


2.865 


-0.38 


-30.81 






57 


36 


41 


318 


4920.5 


2.832 


0.06 


-30.83 






81 


58 


64 


318 


4957.2 


2.851 


-0.42 


-30.84 






56 


34 


39 


318 


4957.5 


2.808 


0.19 


-30.85 






95 


67 


75 


342 


6229.2 


2.845 


-2.94 


-31.47 


19 


10 








342 


6270.2 


2.858 


-2.57 


-31.47 


28 


20 








383 


5068.7 


2.940 


-1.11 


-30.82 






18 


11 


14 


383 


5139.4 


2.940 


-0.59 


-30.83 






44 


27 


21 


383 


5191.4 


3.038 


-0.61 


-30.80 






36 


19 


25 


383 


5232.9 


2.940 


-0.13 


-30.86 






67 


45 


50 


383 


5266.5 


2.998 


-0.43 


-30.83 






47 


29 


34 


383 


5281.7 


3.038 


-0.91 


-30.82 






22 


13 


14 


384 


4787.8 


2.998 


-2.62 


-30.70 


18 










553 


5217.3 


3.211 


-1.03 


-30.71 


88 


95 


10 


6 


7 


553 


5253.4 


3.283 


-1.57 


-31.31 


56 


44 








553 


5324.1 


3.211 


-0.05 


-30.75 






50 


32 


36 


553 


5339.9 


3.266 


-0.62 


-30.71 






21 


13 


13 


553 


5393.1 


3.241 


-0.65 


-30.75 






20 


11 


13 


554 


4736.7 


3.211 


-0.80 


-30.48 


110 




23 


13 


14 


686 


5586.7 


3.368 


-0.07 


-30.71 






42 


27 


28 


686 


5569.6 


3.417 


-0.54 


-30.66 






21 


12 


13 


686 


5572.8 


3.396 


-0.28 


-30.68 






32 


19 


21 


686 


5615.6 


3.332 


0.05 


-30.74 






52 


33 


37 


686 


5624.5 


3.417 


-0.75 


-30.68 


100 




20 


8 


9 


816 


6232.6 


3.654 


-1.21 


-30.67 


61 


44 








816 


6246.3 


3.602 


-0.71 


-30.71 


89 


93 






6 


816 


6400.0 


3.602 


-0.25 


-30.74 






25 


15 


17 


816 


6411.6 


3.654 


-0.62 


-30.71 


98 




13 


8 


9 


984 


4985.2 


3.928 


-0.71 


-30.08 


86 


74 








1031 


5491.8 


4.186 


-2.22 


-30.53 


6 










1062 


5525.5 


4.230 


-1.28 


-30.05 


36 










1087 


5662.5 


4.178 


-0.60 


-30.07 


75 


56 








1087 


5705.4 


4.301 


-1.49 


-30.05 


21 


9 








1092 


5133.6 


4.178 


0.20 


-30.40 






25 


15 


16 


1094 


5074.7 


4.220 


-0.10 


-30.32 


104 


89 


14 


9 




1146 


5364.8 


4.445 


0.15 


-30.21 


108 


94 


16 


11 


9 


1146 


5367.4 


4.415 


0.22 


-30.25 






21 


14 


11 


1146 


5369.9 


4.371 


0.34 


-30.31 






25 


17 


15 


1146 


5383.3 


4.312 


0.42 


-30.38 






30 


21 


20 


1146 


5424.0 


4.320 


0.52 


-30.39 






34 


22 


22 


1164 


5410.9 


4.473 


0.14 


-30.20 




102 


19 


11 


13 


1164 


5415.1 


4.386 


0.41 


-30.31 






29 


18 


17 


1178 


6024.0 


4.548 


0.03 


-30.40 


93 


71 


10 


6 


5 


1179 


5855.1 


4.607 


-1.56 


-30.26 


12 










1195 


6752.7 


4.638 


-1.27 


-30.05 


20 
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Table 5. Atomic data of Fe n lines used in the analysis of the reference stars. Column 1 gives the multiplet number, column 2 the wavelength 
in A, column 3 the lower excitation energy, column 4 the differential gf value with respect to log £(Fe) Q = 7.51, column 5 the employed 
logC 6 value. A radial-tangential approximation to macroturbulence (Gray 1977) was used for the analysis of the Sun, a Gaussian otherwise. 
The equivalent widths (in mA) in columns 6- 10 are of low quality, but only enter the analysis in the derivation of the microturbulence. We 
discourage from using these values for abundance analyses 

mult. A E low logg/auf logC 6 Procyon Gmb 1830 HD 19445 HD 84937 HD 140283 



27 


4416.8 


2.77 


-2.64 


-31.78 




28 


14 


14 


1 1 


27 


4233.1 


2.57 


-2.04 


-31.78 




77 


44 


47 


47 


35 


5136.8 


2.83 


-4.44 


-31.78 


20 










35 


5132.6 


2.79 


-4.17 


-31.78 


35 










35 


5100.6 


2.79 


-4.18 


-31.78 


37 










36 


4993.3 


2.79 


-3.73 


-32.18 


49 










37 


4629.3 


2.79 


-2.38 


-32.18 






20 


19 


17 


37 


4582.8 


2.83 


-3.20 


-32.18 


86 










37 


4555.8 


2.82 


-2.48 


-32.18 




33 


18 


19 


17 


37 


4515.3 


2.83 


-2.53 


-31.88 




32 


14 


15 


14 


37 


4491.4 


2.84 


-2.75 


-31.78 


108 


21 


9 


9 




38 


4620.5 


2.82 


-3.33 


-31.78 


78 










38 


4576.3 


2.83 


-2.95 


-31.78 


96 




6 


6 




38 


4508.2 


2.84 


-2.39 


-31.78 






18 


20 


16 


40 


6516.0 


2.88 


-3.29 


-32.11 


79 










40 


6432.6 


2.88 


-3.61 


-32.11 


67 










40 


6369.4 


2.88 


-4.19 


-32.11 


33 










41 


5284.1 


2.88 


-3.11 


-32.11 


88 


1 1 








42 


5169.0 


2.88 


-1.28 


-32.01 






83 


80 


70 


42 


5018.4 


2.88 


-1.29 


-32.11 




96 


70 


72 


70 


42 


4923.9 


2.88 


-1.53 


-31.91 




87 


61 


65 


58 


43 


4656.9 


2.88 


-3.62 


-32.11 


58 










4(S 


(S084 1 

UUOt. 1 


3 1Q 


-3 84 


-39 1Q 












46 


5991.3 


3.14 


-3.56 


-32.19 


51 










48 


5362.8 


3.19 


-2.57 


-32.19 






8 


8 


7 


48 


5264.8 


3.22 


-3.08 


-32.19 


73 


6 








49 


5425.2 


3.19 


-3.27 


-32.19 


65 










49 


5325.5 


3.21 


-3.21 


-32.19 


68 










49 


5316.6 


3.14 


-1.91 


-31.89 




44 


27 


29 


26 


49 


5234.6 


3.21 


-2.21 


-31.89 




28 


15 


16 


14 


49 


5197.5 


3.22 


-2.27 


-31.89 




23 


13 


13 


11 


74 


6456.3 


3.89 


-2.09 


-32.18 


100 


11 


6 


7 




74 


6416.9 


3.87 


-2.67 


-32.18 


63 










74 


6247.5 


3.87 


-2.33 


-32.18 


89 


6 








74 


6239.9 


3.87 


-3.49 


-32.18 


29 










74 


6149.2 


3.87 


-2.76 


-32.18 


64 











tradiction lies in the choice of Ce values: The larger C(> values 
as proposed by Anstee & O'Mara (1991, 1995) enter our anal- 
ysis via the evaluation of the strong solar Fe i lines where a cer- 
tain exchange between /values and damping constants is pos- 
sible (for more details on our exact choice of damping param- 
eters see Paper II). Thus our set of differential /values differs 
systematically from that employed in the studies by Fuhrmann 
(1998a) and Fuhrmann (1998b, 2003) which employ C 6 values 
calculated using the Unsold (1968) approximation. 

We have no physical explanation for our empirical finding 
that hydrogen collisions are as efficient as Sh = 3. Laboratory 
measurements for the Na i resonance line at H beam energies 
above 10 eV indicate efficiencies as low as Sh~0.01 (Fleck et 
al. 1991). We are convinced that Sh = 3 is not an artefact of 



our modelling, as the same code results in much lower effi- 
ciencies for simpler atoms like Mg or Al (Baumiiller & Gehren 
1996, Zhao et al. 1998). But whereas Mgi and Ah are essen- 
tially shaped by a few strong photoionization cross-sections, 
nearly all Fe i terms are depopulated by photoionization. This 
gives Fe i a unique atomic structure. In spite of the laboratory 
results mentioned above, it can thus not be ruled out that hy- 
drogen collisions among Fe i terms are in fact as efficient as 
the Drawin (1968, 1969) formula predicts (or three times as 
efficient). We note that Sh would even be higher, if we hadn't 
enforced upper-level thermalization above 7.3 eV to account for 
term incompleteness (for details see Paper I). 

Holweger (1996) uses the red giant Pollux (= y Gem) to 
calibrate S H to a value of 0. 1 . This result is not in direct con- 
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Table 6. Different sets of stellar parameters for Procyon. The radius is 
determined from T e g and Mm and has to be compared with the funda- 
mental value of R = (2.07 ± 0.02) Rq. A mP = 100(rf spec - d mP )/d mP 

model r cff logg [Fen/H] diff mass radius A H ip 
[K] ' [M ] [Re] [%] 



1 T cff , logg, [Fe/H]: Allende Prieto et al. (2002) 

2 7" eff : Balmer profiles/log g: astrometry/[Fe/H]: from Fen 

3 r cff : Balmer profiles/log g: Fe i/Fe n (Sh = 3) 

4 r e ff : free parameter, log g: astrometry and Fe i/Fe n (Sh = 3) 



flict with ours as Holweger did not have access to the Bautista 
(1997) cross-sections for photoionization. 

From the point of view of hydrodynamical model atmo- 
spheres it could be argued that the model temperatures in the 
outer layers of static atmospheres are overestimated (Asplund 
et al. 1999). Lower temperatures would predominantly make 
low-excitation lines of neutral species come out stronger, 
thus alleviating the need for strongly thermalizing collisions 
(high values of Sh) to counterbalance photoionization. Using 
2D radiation-hydrodynamical simulations, Steffen & Holweger 
2002 derive granulation corrections for lines arising from low- 
lying levels of neutral species reaching up to -0.3 dex in the 
solar case. More work on combining hydrodynamical with non- 
LTE computations is clearly needed (cf. Shchukina & Trujillo 
Bueno (2001) for first attempts with respect to Fe). 

4.4. Procyon 

Apart from the Sun, a CMi is one of the very few stars for 
which the whole set of fundamental stellar parameters except 
metallicity is known from direct measurements. Allende Prieto 
et al. (2002) review recent determinations of these quantities. 
Its proximity (d-(3.5 + 0.01) pc) and the presence of a white 
dwarf companion (Procyon B) allow to constrain Procyon's 
mass (M = (1.42 ± 0.06)M o ). From the radius (R = (2.01 + 
0.02) R ) and measurements of the integrated flux the effective 
temperature is known to within ~ 100K and probably lies be- 
tween 6500 K and 6600 K. Mass and radius can be combined to 
yield the logarithmic surface gravity, logg = 3.96 + 0.02. The 
metallicity is close to solar. Allende Prieto et al.'s final choice 
for their comparative ID static and 3D hydrodynamical mod- 
elling is T e g — 65 12 K, log g — 3.96 and solar composition. We 
adopt the same value for the surface gravity and derive the ef- 
fective temperature from the Balmer lines which yield an av- 
erage value of r e ff = (6508 + 60) K. The evolutionary tracks 
of VandenBerg et al. (2000) point towards a somewhat higher 
mass of 1.52M , in better agreement with the mass estimate of 
Girard et al. (2000) who also consider ground-based measure- 
ments for the angular separation of Procyon A and B. To be 
able to compare our results directly to those of Allende Prieto 
et al. (2002), we stick to the choice of parameters given above. 



l.oo - 




6540 6560 6580 6600 

wavelength [A] 

Fig. 3. Profile fits to Ha in the FOCES spectrum of Procyon. The re- 
vised temperature (r cff = 6600 K) is in disagreement with the observa- 
tions 



We note, however, that this choice of mass produces a 3 % off- 
set with respect to Hipparcos. 

It turns out that our model is not able to fulfil the ionization 
equilibrium of iron with Sh = 3. With differential non-LTE ef- 
fects in abundance of merely 0.01 dex the unbalance amounts 
to A(Fen-Fei) = 0.07 dex. To remove this discrepancy this 
model would require a gravity of logg = 3.81, many standard 
deviations away from the astrometric result. 
How would one have to tune Sh to remove A(Fen-Fei)? 
Surprisingly, no choice of Sh is able to fulfil the ionization 
equilibrium. Using the largest non-LTE corrections possible 
(Sh = 0) the discrepancy is merely reduced by 0.02 dex. In this 
sense Procyon only documents our failure to use this star as a 
calibrator for Sh- 

Seemingly, Allende Prieto et al. (2002) are more suc- 
cessful with their ID LTE modelling of Procyon's iron spec- 
trum. Their A(Fen-Fei) amounts to 0.02 dex in log e(Fe) 
(log e(Fei)= 7.30 + 0.11 vs. log e(Feii) =7.32 ± 0.08) or 
0.04 dex in \ogg. Upon closer inspection one notices that these 
values are non-differential abundances. Formulated relative to 
the Sun where a discrepancy of 0.06 dex was obtained by 
these authors, the discrepancy in Procyon rises to 0.08 dex 
which is in complete numerical agreement with our LTE result. 
Irrespective of whether one prefers differential analyses or not, 
we conclude that Allende Prieto et al. (2002) do not master iron 
in this star, either (their 3D hydrodynamical modelling results 
in a residual trend of log e(Fe i) with line strength which could 
be removed by including a microturbulence of £ » 0.3 km/s). 

Steffen (1985) investigated different ionization equilib- 
ria in Procyon and drew the conclusion that the fundamen- 
tal temperature needs to be revised upwards by more than 
200 K. Alternatively, the gravity would have to be lowered to 
log g =i 3.55 which can be ruled out by astrometry. In the frame- 
work of non-LTE modelling, a similar temperature correction 
would be required to bring [0 1 6300] into concordance with the 



1 


6512 


3.96 


-0.09 


1.42 


2.13 


-3.1 


2 


6508 


3.96 


-0.03 


1.42 


2.13 


-3.0 


3 


6508 


3.81 


-0.10 


1.42 


2.13 


+ 15.0 


4 


6600 


3.96 


-0.03 


1.42 


2.06 


-0.1 
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Table 7. Stellar parameters derived using Fei/n in non-LTE as a gravity indicator with Sh = 3. The spectroscopic distance rf spec is calculated 
from log ;r S pe C = 0.5 ([g] - [M]) - 2 [T sB ] - 0.2 (V + BC + A v + 0.25), where [X] = \og(X/X e ). A H ip = 100 (d spec - d mP )/d mP . The oxygen 
abundance is derived by means of profile analysis of the IR triplet in non-LTE (Reetz 1999), the magnesium abundance from profiles of weak 
optical lines in non-LTE (AA 4571, 4702, 4730, 5528, 5711 A, Zhao et al. 1998) 



object 


T e ff 


log 8 


[Fe/H] 


[O/Fe] 


[Mg/Fe] 


mass 


A v 


BC 




dm? 


Amp 




[K] 




NLTE 


NLTE 


NLTE 


[Mo] 


[mag] 


[mag] 


[pc] 


[pc] 


[%] 


HD 103095 


5070 


4.66 


-1.36 


+0.63 


+0.21 


0.62 


0.00 


-0.32 


8.86 


9.16 


-3.2 


HD 19445 


6032 


4.40 


-2.08 


+0.68 


+0.49 


0.67 
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IR triplet lines (Reetz 1999). We can follow along this track 
and seek the temperature at which the Sh = 3 model brings 
loge(Fei) into concordance with loge(Fen) without contra- 
dicting the astrometric distance constraint. This temperature is 
found to be T e g = 6600 K or 90 K higher than our best spectro- 
scopic estimate based on Balmer profiles. As can be appreci- 
ated from inspecting Fig. 3 this temperature is in conflict with 
the Balmer profile analysis. It is, however, in not in disagree- 
ment with the gravity inferred from the Mg ib lines when em- 
ploying the non-LTE model atom for magnesium presented by 
Zhao et al. (1998). This is because the strong-line method is 
not very sensitive to changes in temperature. Thus, while dis- 
quieting on the whole, our result is an improvement over the 
analysis of Steffen (1985). Higher temperatures for Procyon 
certainly seem possible in view of the Stromgren photometry 
of Edvardsson et al. (1993, r eff = 6705K) and the IRFM tem- 
perature of Alonso et al. (1996, = 6631 K). From the new 
self-broadening theory for Balmer lines presented by Barklem 
et al. (2000) higher temperatures (T e g = 6570 K) are possible, 
if the zero-point offset for the Sun is properly accounted for. 

5. Stellar parameters from the S H = 3 model 

Having calibrated S H , we reanalyse the iron spectra of the ref- 
erence stars with our non-LTE model for iron. The final para- 
meters are given in Table 7. 

Fig. 4 compares the stellar parameters in non-LTE with 
their LTE counterparts. In particular, the scatter is reduced 
in going from LTE to non-LTE, as differential non-LTE ef- 
fects for HD 103095 bear the opposite sign to those of all the 
other stars. Non-LTE corrections are largest for HD 84397 and 
HD 140283. A positive slope of a linear regression line in this 
diagram would thus indicate underestimated non-LTE correc- 
tion and vice versa. This effect will be seen in the analyses of 
Thevenin & Idiart (1999) and Fulbright (2000) below. 
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Fig. 4. A H ip as a function of the astrometric distance, both in non-LTE 
(top) and LTE (bottom). The grey wedge indicates the typical uncer- 
tainty of the Hipparcos parallaxes. When non-LTE gravities are con- 
sidered the star-to-star scatter is significantly reduced 



valid comparison, we use the published stellar parameters, but 
rederive bolometric corrections and masses. Since an estimate 
for the [a/Fe] ratio is needed for interpolating among appropri- 
ate evolutionary tracks, we only consider studies which supply 
an a-element abundance alongside [Fe/H]. We therefore disre- 
gard e.g. the study of Soubiran et al. (1998). 

Particular attention is paid to studies which also undertook 
kinetic equilibrium calculations for iron. As can be appreciated 
from comparing the upper panel of Fig. 4 to Fig. 5, all studies 
using Fe i/Fe n are an order of magnitude less precise than our 
calibration. 



5.1. Comparison with other studies 

Based on our newly established methodology, we can now 
compare the resulting spectroscopic distance scale with that 
of other studies. We compare the residual offset with respect 
to the Hipparcos astrometry and the star-to-star scatter for the 
- admittedly - small sample of reference stars. For a fair and 



5.1.1. Gratton et al. (1999) 

These authors employ a rather simplistic model atom of Takeda 
(1991), treat the Fei photoionization and background opaci- 
ties empirically and calibrate Sh by means of RR Lyrae stars. 
As Sh is determined to be 30, thermal populations result for 
their stars which is why the LTE stellar parameters of Gratton 
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et al. (1996) are still the basis for the analysis of Carretta et 
al. (2000). Due to the dominant influence of oxygen on the 
morphology of the evolutionary tracks, we use their published 
[0/Fe] NLTE as [ar/Fe]. 

The topmost panel of Fig. 5 shows the resulting spectro- 
scopic distance scale in comparison with the astrometric dis- 
tances of Hipparcos. Since the temperatures these authors de- 
rive via the IRFM are in excellent agreement with ours, it is 
surprising to see that the gravities are somewhat overestimated 
rather than underestimated. From the work of Fuhrmann et al. 
(1997) we would have expected the LTE gravities to be sig- 
nificantly too low (cf. Introduction). Since we have shown this 
effect to be mainly due to an inappropriate choice of damping 
parameters Ce, we assume that a similar effect is at work in 
the analyses of Gratton et al. (1996). More worrying than the 
12% offset is the large scatter, e.g. in the multiple analysis of 
HD 140283 which amounts to nearly 30 %. 

5.1.2. Thevenin & Idiart (1999) 

In an attempt to model the kinetic equilibrium of iron ab initio, 
Thevenin & Idiart (1999) put together a comprehensive and up- 
to-date model atom which is quite comparable to the one used 
by us. Photoionization is treated according to the calculations 
of Bautista (1997), UV line opacities are crudely represented 
by a scaling rule for continuous opacities. In combination with 
their neglect of hydrogen collisions (Sh = 0), these assumptions 
lead to very large non-LTE corrections. For example, the grav- 
ity correction for HD 140283 is calculated to be 0.54 dex. We 
use the larger value of [Mg/Fe]NLTE and [Ca/Fe]NLTE fr° m Idiart 
& Thevenin (2000) as an estimate of [a/Fe] . 

On the temperature scale employed by these authors, the 
non-LTE corrections are clearly overestimated, severely so in 
the case of HD 84937 (A H ip = -34 %). While the scatter is com- 
parable to that of Gratton et al. (1999), the offset is significantly 
larger: the sample distances are incompatible with Hipparcos at 
the 2cr level. 

These authors also analyse Procyon. At a temperature of 
T e s = 663 1 K the ionization equilibrium of iron is reached at a 
gravity of log g = 4.00 without departures from LTE. Judging 
from our own calculations (cf. Sect. 4.4) this modelling suc- 
cess is not surprising and basically the result of choosing 
a high effective temperature. The derived parameters yield 
A H ip = -3.3%. Including this star into the sample, the mean 
offset would be reduced to 17 % with the scatter increasing to 
12%. 



5.1.3. Fulbright (2000) 

This analysis of 168 mostly metal-poor stars rests almost en- 
tirely on iron as a plasmadiagnostic tool: the temperature de- 
termination is based on the LTE excitation equilibrium of 
Fei, gravity and metallicity are inferred from the LTE ion- 
ization equilibrium of Fei/Fen. These steps are iterated un- 
til log e(Fe i) matches log e(Fe n) to "within 0.03 dex in most 
cases". At best, we can therefore expect an accuracy for log g 
of 0.05 dex. We use [Mg/Fe] LTE as a proxy for [a/Fe]. The 
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Fig. 5. Same as Fig. 4, here for the studies discussed in the text. Note 
the different y-axis scales. Except for the work of Fuhrmann (1998b, 
2000), all studies show large star-to-star scatter and/or offsets with 
respect to Hipparcos 



masses we derive from Fulbright's stellar parameters show 
a wide range: HD 103095 would be a very low-mass ob- 
ject (M = 0.45 M Q ), HD 140283 a star of nearly solar mass 
(M = 0.95 M ). This might indicate unrealistic combinations of 
r eff and logg. 

While the sample mean comes out relatively close to 
Hipparcos, the star-to-star scatter is very large. The offset 
for HD 140283 is 40%. Statistically speaking, the gravities 
are only know to within 0.25 dex. Just like in the case of 
Thevenin & Idiart (1999) not a single star falls in the lcr un- 
certainty range of the Hipparcos astrometry. All this means that 
Fulbright's methods fail to give consistent results for metal- 
poor F and G stars. 



5.1.4. Fuhrmann (1998b, 2000) 

These studies are methodologically closest to our own, in the 
sense that Fuhrmann employs the same model atmosphere 
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and the same Balmer profile temperature scale. The pressure- 
broadened wings of Mgib in LTE are used as a gravity in- 
dicator with the remarkable result of A H i P = (2 + 5) % for 100 
stars. Balmer profiles and Mgib lines are practically orthogo- 
nal methods to determine r eff and log g which helps to prevent 
disadvantageous propagation of errors. Masses are interpolated 
among tracks of Bernkopf (1998). The stellar parameters of 
HD 140283 are revised with respect to Fuhrmann (1998a) and 
are private communication. 

When using the tracks of VandenBerg et al. (2000) the re- 
sulting distances are slightly shorter than with the Bernkopf 
(1998) tracks. The internal consistency is, however, still very 
impressive. Fuhrmann's work clearly demonstrates the capa- 
bilities of a careful, differential analysis in LTE. 

We note in passing that the revision of the stellar mass for 
Procyon (down to M = 1.42M ) worsens the spectroscopic re- 
sult published by Fuhrmann (1998b) to A raP = -9.6%. Thus 
Procyon still is a paramount test case for any spectroscopic 
analysis. 

6. Extreme applications: HE 0107-5240 & 
CS 29497-004 

In the meantime, our model has been applied to a cou- 
ple of new extreme halo objects: HE 0107-5240, the most 
iron-deficient star currently known (Christlieb et al. 2002) 
and CS 29497-004, a new highly r-process enhanced star 
(Christlieb et al. 2003). Since both objects are on the upper red 
giant branch (RGB), non-LTE corrections to gravity are larger 
than in the case of the reference stars discussed above. 
Both stars have effective temperatures around 5100 K and sur- 
face gravities around 2.2. Despite their vastly different metal- 
licities ([Fe/H] = -5.3 vs. -2.6), non-LTE corrections to grav- 
ity are very similar and amount to +0.3 dex. They are similar 
because two competing processes cancel each other. On the one 
hand, the Fe i departure coefficients of HE 0107-5240 are more 
extreme than those of CS 29497-004; on the other, the Fe i 
lines in HE 0107-5240 are much weaker and therefore origi- 
nate deeper in the atmosphere. As test calculations have shown, 
hydrogen collisions are still important to consider at these RGB 
gravities. 

Further calculations indicate that non-LTE effects in Fei 
reach a plateau value below a metallicity comparable to that 
of the most metal-poor globular clusters found in our Galaxy. 
Globular cluster giants are thus ideal targets to verify the pre- 
dictions of our calculations. It is encouraging to see that at 
intermediate metallicities these predictions are in good agree- 
ment with the mild Fe i/Fe n imbalance found in the LTE ana- 
lyses of giants in M5 and M71 by Ramirez & Cohen (2003). 

7. Conclusions 

We have presented a procedure to determine the poorly-known 
efficiency of collisions with neutral hydrogen in our kinetic 
equilibrium calculations for the formation of Fe i and Fe n lines 
in the atmospheres of solar-type stars. This procedure rests 
on exploiting the astrometric constraint of Hipparcos paral- 
laxes. The scaling factor Sh is determined to be 3 which makes 



hydrogen collisions the most important atomic process coun- 
terbalancing the otherwise overwhelming photoionization in 
metal-poor stars. We emphasize, however, that our goal has not 
been to determine the absolute strength of hydrogen collisions. 
Rather, our analysis is aimed at removing biases that are present 
in the analysis of metal-poor stars when using the iron ioniza- 
tion equilibrium in LTE using ID model atmospheres. 

For a representative sample of local halo stars the combina- 
tion of Balmer profile temperatures and non-LTE iron ioniza- 
tion equilibrium gravities succeeds in fulfilling the astromet- 
ric constraint. Other authors (Gratton et al. 1999, Thevenin & 
Idiart 1999, Fulbright 2000) present results an order of magni- 
tude less accurate. This underlines our claim that hydrogen col- 
lisions must be considered in kinetic equilibrium calculations. 
The independent method to derive \ogg used by Fuhrmann 
(1998b, 2000) yields comparable results. 

Our method cannot establish the ionization equilibrium 
of Procyon at the fundamental parameters r e ff = 6510K and 
logg = 3.96. An upward revision of Procyon's temperature by 
90 K ( 1 . 3 %) to r eff = 6600 K is indicated instead, in agreement 
with the radius determination based on its bolometric magni- 
tude. This temperature is, however, in disagreement with the 
spectroscopic constraint of the Balmer profiles. 

Using our carefully calibrated model of Fe i/Fe n, cool stars 
over the whole range of metallicities encountered in the Galaxy 
can now be analysed in a homogeneous way with reduced 
methodological biases without resorting to trigonometric par- 
allaxes. 
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